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1 Introduction 



The search for the standard model (SM) HJ-SlI Higgs boson is one of the main goals of the Large 
Hadron Collider (LHC) physics programme. The discovery of the SM Higgs boson would 
shed light on the spontaneous electro weak symmetry breaking mechanism [$-^|. To date, ex- 
perimental searches for this particle have yielded null results. Limits at 95% confidence level 
(CL) on its mass have been placed by experiments at LEP, wh > 114.4GeV IITOt , the Teva- 
tron, niH ^ (162-166) GeV CIl, and ATLAS, nin i (145-206), (214-224), (340-450) GeV IlMll. 
Precision electroweak measurements, not taking into account the results from direct searches, 
indirectly constrain the SM Higgs boson mass to be less than 158 GeV [15J. 

This letter presents a search for the SM Higgs boson in the decay mode H — )■ ZZ — )■ £+£^t+t^, 
where I is either }i or e. One Z is required to decay either into or e+e^, and the second 

Z into T+T^ in four possible final states: ThT^, t^t^, TeT^, and TeT^„ where Th represents a r 
decaying hadronically, and T^,, and Te indicate taus decaying into muons and electrons respec- 
tively. The present measurement complements the search in the H — > ZZ — ?> 4£ channel fl^. 
The presence of four leptons in the final state provides a clean signature with only a small 
contribution from background processes. The major irreducible background contribution is 
non-resonant ZZ production. The most important reducible background contributions are Z 
and WZ production in association with jets, and tt production. The final states tt — >• t^t^;, TgTe 
are not considered, as they are accounted for in the H — > ZZ — > 4£ Higgs search |1T6| . The cross 
sections for the Higgs boson production mechanisms and decay branching fractions, together 
with their uncertainties, are taken from Ref. I|T7| and are derived from Refs. Iil8l - I37| . 

The analysis is based on data from proton-proton collisions at v^s = 7 TeV, corresponding to an 
integrated luminosity of 4.7 fb^^ collected with the Compact Muon Solenoid (CMS) detector at 
the LHC in 2011. This is the first Higgs boson search performed in the H — > ZZ — )■ I^t^t^ 
channel. 



A detailed description of the CMS detector can be found elsewhere [38J. The central feature 
of the CMS apparatus is a superconducting solenoid, of 6 m internal diameter, providing a 
magnetic field of 3.8 T. Within the field volume are the silicon pixel and strip tracker, the crystal 
electromagnetic calorimeter (ECAL), and the brass /scintillator hadron calorimeter. Muons are 
measured in gas-ionization detectors embedded in the steel return yoke. 

CMS uses a right-handed coordinate system, with the origin at the nominal interaction point, 
the :t-axis pointing to the centre of the LHC ring, the y-axis pointing up (perpendicular to the 
LHC plane), and the z-axis along the counterclockwise-beam direction. The polar angle, Q, is 
measured from the positive z-axis and the azimuthal angle, ^, is measured in the x-y plane. 
Variables used in this analysis are the pseudorapidity, f] = — In [tan (0/2)], and the transverse 



The ECAL is designed to have both excellent energy resolution and high granularity, prop- 
erties that are crucial for reconstructing electrons and photons produced in T-lepton decays. 
The ECAL is constructed with projective lead tungstate crystals that provide coverage in pseu- 
dorapidity \ri\ < 1.479 in a barrel region and 1.479 < \t]\ < 3.0 in two endcap regions (EE). A 
preshower detector consisting of two planes of silicon sensors interleaved with a total of 3 Xq of 
lead is located in front of the EE. The energy resolution is 3% or better for the range of electron 
energies relevant for this analysis. 



2 CMS Detector 



momentum, pj 
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3 Event Selection and Monte Carlo Samples 



The inner tracker measures charged particle tracks within the range \r]\ < 2.5. It consists of 
1 440 siHcon pixel and 15 148 silicon strip detector modules, and provides an impact parameter 
resolution of ~15f/m and a transverse momentum resolution of about 1.5% for 100 GeV par- 
ticles. The reconstructed tracks are used to measure the location of interaction vertices. The 
spatial resolution of the reconstruction is ~25 ^m for vertices with more than 30 associated 
tracks Il39l . 

The muon barrel region is covered by drift tubes, and the endcap regions by cathode strip 
chambers. In both regions, resistive plate chambers provide additional coordinate and timing 
information. Muons are reconstructed in the range \ri\ < 2.4, with a typical pj resolution of 
f=:il%for Pt ~ 40 GeV. 

3 Event Selection and Monte Carlo Samples 

At the trigger level, the selected events are required to have at least two muons, one with pj > 
13 GeV (pt > 17 GeV for the end of the data-taking period when the instantaneous luminosity 
was highest) and the other with pj > 8 GeV, or at least two electrons, one with pj > 17 GeV 
and the other with pj > 8 GeV. 

Algorithms for identifying muons and electrons, collectively referred to as leptons, are based on 
the tracker, the muon systems and the calorimeters f40', 14111. Since the ZZ final state is expected 
to have only a small contribution from background processes, the algorithms are tuned to max- 
imize the lepton-reconstruction efficiency, resulting in an increased lepton-misidentification 
rate. A particle flow (PF) technique [42] is used to form lepton-isolation quantities and is 
also used for Th reconstruction. In the PF approach, information from all subdetectors is com- 
bined to reconstruct and identify particles produced in the collision. The particles are classified 
into mutually exclusive categories: charged hadrons, photons, neutral hadrons, muons, and 
electrons. These particles are used to reconstruct Th with the "hadron plus strip" (HPS) algo- 
rithm |43J that is designed to optimize the performance of identification and reconstruction 
by considering specific Th decay modes. The neutrinos produced in all t decays escape detec- 
tion and are ignored in the Th reconstruction. The algorithm provides high Th identification effi- 
ciency, approximately 50% for the range of Th energies relevant for this analysis, while keeping 
the misidentification rate for jets at the level of ~1%, that is factor of three to four times lower 
with respect to other available algorithms [|44|. 

Events are required to have at least one Z ^ candidate, denoted by Zi, with the leptons 
of opposite charge, one with pj > 20 GeV and another with pj > 10 GeV, and with \tj\ < 2.4 
for the muons and |?/ 1 < 2.5 for the electrons. Both leptons are required to have a combined PF 
relative isolation J^^ < 0.25, which is defined as: 

jPF ^ ^pCharged ^ ^^^^^y ^ ^neutral _ 3 5 ^ p^^ ,0)) / p'j, (1) 

where p*^"^^*^*^ is the scalar sum of the charged hadrons px, and Ej and £3^'^"'''^l correspond, 
respectively, to the sum of the transverse energies of the photons and neutral hadrons, all mea- 
sured in the isolation cone of AR < 0.4 around the lepton direction, where AR = ^ (A//)^ + (A^)^. 
The contribution from neutrals is corrected for the effect of overlapping pp interactions. The 
value of the correction is estimated by scaling the sum of the pj of all charged particles not 
associated with the vertex in the isolation cone (pj^) by a factor 0.5, which is approximately 
the ratio of neutral to charged hadron yields in hadronization processes. 

For the second Z, denoted by Z2, the selection requirements depend on the final state. If the 
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final state is T^^Te, the lepton pj values are required to exceed lOGeV. The remaining criteria 
are identical to those for Zi. Since t^s have much larger misidentification rates than the other 
leptons, the isolation requirement based on for the muons and electrons in the final states 
T^Th and TeTh is changed to 0.15 and 0.1, respectively. In a study of inclusive Z — > tt produc- 
tion 1451 , it was demonstrated that modifying the muon and electron isolation requirements is 
a more effective way to reduce background in such final states than requiring tighter isolation 
on Th. The are required to have pj > 20 GeV and |?/ 1 < 2.3, and to pass a loose HPS working- 
point requirement. If the Z2 decays to both are required to pass a medium working 
point of the HPS algorithm. The loose (medium) working point requires the scalar sum over 
the charged hadrons px arid the neutral hadrons Ej in the isolation cone, to be less than 2 GeV 
(1 GeV). The loose (medium) working point corresponds to a probability of approximately 1% 
(0.5%) for jets to be misidentified as Th. Using the medium instead of loose working point leads 
to a decrease in the t^, reconstruction efficiency from ?»50% to ?»40%. 

The visible invariant mass of the reconstructed Z2 — > t^t^ is required to be 30 < nirr < 
80 GeV, and this criterion is used for most of the final states. The upper bound reduces contri- 
butions from Z2 — ?► ££, where a muon or an electron is not well reconstructed, and misidentified 
as a Th. For the Z2 — > TeT„ final state, the upper bound is increased to 90 GeV, as this state is 
not produced in Z2 — > ££ decays. Leptons from the same Z are required to be separated by 
AR > 0.4 for Zi, and by AR > 0.5 for Z2. The two reconstructed Zi and Z2 are required to be 
separated by AR > 0.5. 

A set of Monte Carlo (MC) event samples is used to simulate signal and background events. 
The Drell-Yan background, £^i^ in association with jets, is simulated with the next-to-leading 
order (NLO) MC generator powheg 2.0 ll46H48ll . The QCD multijet, W and diboson WZ back- 
grounds are simulated with pythia 6.424 l|49| . The ZZ background is simulated with pythia 
6.424 and MadGraph [50J. The tf samples are simulated with MadGraph. The r-lepton de- 
cays are generated with TAUOLA [51]. The Higgs boson samples are generated using POWHEG 
2.0, which incorporates NLO gluon fusion (gg — ?> H) and vector-boson fusion (qq — > qqH). All 
events are processed through a detailed simulation of the CMS detector based on GEANT4 1152 1 
and reconstructed with the same algorithms that are used for data. 



4 Background Estimates and Systematic Uncertainties 

The major irreducible source of background to the H — > ZZ — > £+£^t+t^ process is from 
SM ZZ — > i^£^T^T^ production. The ZZ contribution is estimated from data by scaling the 
prediction from simulation to the well measured inclusive Z production cross section. The 
number of estimated ZZ events, N^^, can be written as: 

n-SM A 

Ts/est ^ ^obs .^.flZZ ^2) 
JNfzZ ^^Z y^^' 

where N*^ is the number of observed events from inclusive Z production, Az is their esti- 
mated acceptance from a MC simulation, including all selection requirements, and rescaled by 
measured data/MC correction factors, Azz is the acceptance for ZZ events, is the SM cross 
section for inclusive Z production, and is the SM cross section for ZZ production calculated 
withMCFM [53 J. 

The other major background contributions arise from the production of Z and WZ in associa- 
tion with jets, as well as tt and QCD multijet production. The latter two backgrounds are small. 
In all these cases, a jet or non-isolated lepton is misidentified as a Th, Te or t^. The probability 
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5 Results 



for jets to be misidentified as Th is measured using i^i^ZhTh events in data in which the Zi 
passes all selection requirements, but no requirement is applied on isolation, and the two Th 
candidates are required to have the same charge. This region is dominated by Z+jets events. 
The Th misidentification rate is defined as the ratio of the number of Th candidates that pass 
the HPS loose or medium working-point requirements, to the initial number of Th candidates, 
and is measured as a function of the pj for each Th. To estimate the number of background 
events in the signal region, the measured misidentification rate is applied to events that pass 
all selection requirements, including the opposite-charge requirement for the Z2, but requiring 
the Th candidates to not be isolated. 

The misidentification rate for Te and in the jifiT^Ze and eeT^Te final states is estimated using 
events in which the Zi passes all selection requirements, and the event contains an additional 
muon or electron. No isolation requirements are applied to it. The misidentification rate is 
defined for Te and t^, in the same way as described above for Th and applied to i-ijiT^iTe and 
eeZfiTe events that pass all the selection requirements, but requiring Te or t^, to not be isolated. 
Isolated muons and electrons from H ZZ — )■ Ai and ZZ — > Ai production can also be 
misidentified as Th. Events are rejected if they are also identified as ZZ — ?► A£ events with 
criteria described in Ref. |16|. 

Theoretical uncertainties on the Higgs boson cross section (17-20%) and branching ratio (2%) 
are taken from Ref. IITtI . Recent studies [|T7l[54ll55l show that current MC simulations do not 
describe the correct Higgs boson mass line shape above ~300 GeV. This effect amounts to an 
additional uncertainty on the theoretical cross section, and hence on the limits, of about 4% at 
mil = 300 GeV and 10-30% for nin of 400-600 GeV. The main uncertainty on the estimate of 
the ZZ background arises from the theoretical uncertainty on the ZZ production cross section 
and is taken from Ref. ||56 [ . The uncertainties on the other backgrounds, Z+jets, WZ+jets, and 
tt reflect the uncertainties on the measured values of the misidentication rates and the limited 
statistics of the control regions in the data. The uncertainty on integrated luminosity of the 
data sample is 4.5% |57|. Systematic uncertainties on trigger efficiency (1%) and on lepton 
identification efficiency and isolation are evaluated from data. The uncertainties associated 
with lepton identification and isolation are 1-2% for muons and electrons, and 6-7% for Th. 
Uncertainties on energy scales, 3% for Th and 1-2.5% for electrons, contribute to variation in the 
shape of the mass spectrum. 

5 Results 

Ten i^i^T^T^ candidates are observed in eight search channels, while 11.60 ± 0.54 (stat.) ± 
1.62 (syst.) background events are expected. Table [l] compares the estimated number of back- 
ground events to the number of events observed in the signal region. The distribution of the 
reconstructed invariant mass summed over all eight i^l^r^r^ decay channels is shown in 
Fig. [1] The shape of the background is taken from the MC simulation, with each component 
normalized to the corresponding estimated value from Table [l] The expected mass distribu- 
tions for the SM Higgs boson with a mass of = 200 GeV and 400 GeV are also shown in 
Fig. [ij The reconstructed masses are shifted with respect to the generated values by ~ 30% due 
to the undetected neutrinos in t decays. As a result, the H — > ZZ — )■ £^£^t^t^ mass resolution 
is 10-15%, depending on the final state, and is almost independent of wh- 

The product of the acceptance and branching fraction for the individual T-decay channels 
ranges between 0.01-0.02 for mn = 200 GeV, and increases by a factor of three to four for 
mil = 400-450 GeV. This behaviour is expected. The final-state leptons produced in t decays 
of more massive Higgs bosons have higher momenta than those from direct Z ^ i£ production 
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Table 1: The estimated yields of ZZ and other background events obtained from data, as de- 
scribed in the text, are shown for each decay channel and are summed in the total background 
yield ("Total backgr."), and compared to the number of events observed in the signal region. 
The total uncertainty is the sum in quadrature of statistical and systematic uncertainties. The 
number of signal events expected for the SM Higgs boson with a mass of Wh = 200 GeV is also 
shown. 



Decay 
channel 




Other 
backgrounds 


Total 
backgr. 


ma 
200 GeV 


Observed 




0.79 ± 0.09 


0.76 ± 0.31 


1.55 ±0.32 


0.17 





eeThTh 


0.75 ± 0.09 


0.73 ± 0.32 


1.48 ±0.33 


0.15 


1 


eeTeTh 


1.12 ±0.13 


0.99 ± 0.34 


2.11 ±0.36 


0.25 


3 




1.20 ±0.14 


0.31 ± 0.29 


1.51 ±0.32 


0.26 


3 




1.08 ±0.13 


0.67 ±0.36 


1.75 ±0.38 


0.23 


2 




0.94 ±0.10 


0.41 ±0.16 


1.35 ±0.19 


0.20 





eeZeT^ 


0.51 ± 0.06 


0.58 ± 0.42 


1.09 ±0.42 


0.11 







0.58 ± 0.07 


0.18 ±0.18 


0.76 ± 0.22 


0.12 


1 


Total 


6.97 ± 0.84 


4.63 ± 1.49 


11.60 ± 1.71 


1.49 


10 
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Figure 1: The four-lepton reconstructed mass summed for all i^i^T^T^ final states. The data 
corresponds to an integrated luminosity of 4.7 fb^^. Points represent the data, shaded his- 
tograms represent the background and hashed histograms represent the signal expectations 
for two Higgs boson masses. The background shapes are taken from MC simulation and are 
normalized to the values obtained using control data samples, as described in the text. 
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6 Summary 



and the selection requirements and the lepton reconstruction become more efficient at larger 
oth- While the cross section decreases with increasing jmh, the expected number of Higgs bo- 
son events selected by this analysis in 4.7fb^^ of pp collisions at 7TeV is 1.4-1.5 in the range 
200 < niii < 400 GeV, and decreases rapidly at higher masses. 




200 300 400 500 600 

(GeV) 

Figure 2: The expected and observed 95% CL upper limits on the product of the Higgs boson 
production cross section and decay branching fraction normalized to the SM expectation, de- 
noted by cr/ crc,M, as a function of Wh- The one- and two-standard-deviation variations from the 
expected limit are also shown by the green and yellow bands. 

In Fig. |2] the expected and observed upper limits at 95% CL on the product of the Higgs boson 
production cross section and decay branching fraction normalized to the SM expectation are 
presented as a function of mn- The limits are calculated with the modified frequentist construc- 
tion CLg M58l - l60| based on the shape of the £^£^t^t^ invariant mass distributions by including 
all eight individual channels in the likelihood combination. The green and yellow bands rep- 
resent the one- and two-standard-deviation variations from the expected limit. The systematic 
uncertainties are introduced in the form of nuisance parameters with log-normal probability 
density functions. The upper limit on the cross section is approximately a factor four to twelve 
larger than the SM Higgs boson production cross section in the range of 190 < wjh < 600 GeV. 



6 Summary 

A search for the standard model Higgs boson has been performed in the decay mode H — > 
ZZ — )■ £+£^T+T^ using CMS data corresponding to an integrated luminosity of 4.7fb^^. No 
evidence is found for a significant deviation from the background expectation. An upper limit 
four to twelve times larger than the predicted value is set at 95% confidence level for the prod- 
uct of the standard model Higgs boson production cross section and decay branching fraction 
in the mass range 190 < niu < 600 GeV. This is the first Higgs boson search performed in the 
H ^ ZZ ^ £+e-T+T- channel. 
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